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Abstract 



Motivated by the paper by J. B er ger and K. Rubinstein [BeRu 



and 

other recent studies |GiPh|, [LuPal|, [LuPa2|, we analyze the Ginzburg- 
Landau functional in an open bounded set Q. We m ainly discuss the bifur- 
cation problem whose analysis was initiated in | Od and show how some of 
the techniques deve loped by the first author in the case of Abrikosov's su- 
perconductors [Du can be applied in this contex t. In the case of non sim- 
ply connected domains, we come back to [BeRu| and |HH001|, |HH002| 
for giving the analysis of the structure of the nodal sets for the bifurcating 
solutions. 



1 Introduction 



X 



1.1 Our model. 



Following the paper by Berger- Rubinstein | BcRu| ] , we would like to understand 
the minima (or more generally the extrema) of the following Ginzburg-Landau 
functional. In a bounded, connected, regulaij^, open set fl C and, for 
any A > and k > 0, this functional Ga,k is defined, for u S H^{^l;(U) and 
A e Hl^{IR^-IR^) such that rot A e L^, by 



GxAu.A) = (A(-|^/|2 + + |(V ~ iA)u\-'))dx^ ■ dx2 

+K-^A"^ /jj2 IrotA — H(,\^dxi ■ dx2 ■ 



(1.1) 



Here, for A = (^1,^2), rot^ = dxiA2 — 8x2^1, div A = Sx^^i +922^2 and fl"e is 
a function on IR^ (or more generally some function in L^(iR^)). Physically 
He represents the exterior magnetic field. 
Let Ae be a solution of 

rotAe = He „^ 

divAe-0. ^ > 



^with C°° boundary 
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It is easy to verify that such a solution exists by looking for in the form 
Ae = {—dx^tf'e, dxitpe)- We have then to solve Atp^ = and it is known to be 
solvable in S'{IR^) n C^{IR^) (or in S'{IB?) n Hl^{lR^) if e L^ilR^)). Of 
course A^ is not unique but we shall discuss about uniqueness modulo gauge 
transform later and at the end this is mainly the restriction of A^ to which 
will be considered. 

We shall sometimes use the identification between vector fields A and 1-forms 

LOA- 



When analyzing the extrema of the GL-functional, it is natural to first an- 
alyze the corresponding Euler-Lagrange equations. This is a system of two 
equations (with a boundary equation) : 

(GL)i -{V -iAfu + Xu{\u\^ -1) = , in 17 , 

{GL)2 rot*(rot^ - He) = Xk,-^ Im [u • (V - iA)u] ■ In , (1.3) 

{GL)3 (V - iA)u • = , indn . 

Here z/ is a unit exterior normal to Oil. The operator rot* is defined by rot*/ := 
idxj,-dxj). 

Moreover, without loss of generality in our problem, we shall add the condition 

{GL)i divA = in Q. (1.4) 

One can also assume if necessary that the vector potential satisfies 

A-z7=0, (1.5) 

on the boundary of Q, where is a normal unit vector to dfl. 

Let us briefly recall the argument. One would like to find 9 in C°°{Q,) such 



that A = A + dO satisfies (1.4) and (1.5). One can proceed in two steps. The 



first step is to find a gauge transformation such that (1.5) is satisfied. This is 
immediate if the boundary is regular. 
We now assume this condition. 
The second step consists in solving 

AO = -divA in Q, 
1^ = 0, ondn. 

This is a Neumann problem, which is solvable iff the right hand-side is orthogo- 
nal to the first eigenfunction of the Neumann realization of the Laplacian, that 
is the constant function x i-^ 1. We have only to observe that div A dx = 



if (1.5) is satisfied. 



An important remark is that the pair (0, Ae) is a solution of the system. 
This solution is called the normal solution. Of course, any solution of the form 
(0, Ae + V(j)) with (p harmonic is also a solution. 
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Remark 1.1 . 

Note also that the normalization of the functional leads to the property that 

GA,«(0,Ae) =0. (1.6) 
The first proposition is standard. 
Proposition 1.2 . 

Iffl is bounded, the functional G x^k, admits a global minimizer which is a solution 
of the equation. 



We refer to |DGP|, for a proof together with the discussion of the next subsec- 
tion. 

1.2 Comparison with other models 

Let us observe that there is another natural problem which may be considered. 
This is the problem of minimizing, for (u,A) e H^{n,€) X H^{n,IR^), the 
functional G^^ defined by 

G'U^.A) =^(A(-|^|2 + i|«|4) + |(V-zA)«|2) 

+K^A"^ IrotA - i7e|2 • da;2 . ^ ' ' 

This may lead to a different result in the case when U, is not simply connected. 
According to discussion with Akkermans, this is the first problem which is the 
most physical (See also the discussion in the appendix). 

A comparison between Ga.k and G^'^ where Z? is a ball containing and G^^'^ 
is defined by 

Gl^{u,A) =^(A(-H2 + i|^|4) + |(V-a)zi|2) 
+k2a-i |rot a - iJep dxi ■dX2 ■ 

is useful. If b is given with support outside of the ball D, it is easy to see 
(assuming that b is regular) that there exists a with support outside D such 
that rot a = b. It is indeed sufficient to take the usual transversal gauge 



fli = —X2 / sb{sx)ds , a2 = xi / sb{sx)ds . (1-9) 
Jo Jo 

This shows that, for any D containing 17, wc have 

inf Ga,.(w, A) = inf G^'f (u, A) . (1.10) 

In particular it is enough to consider minimizing sequences (u„, + a„) where 

suppa„ C D and D is a ball containing 17. The proof of the existence of 
minimizers is then greatly simplified. 

Finally, it is natural to think that one can replace D by 

n-~nuiO,, (1.11) 

■^This is at least clear whe n Q is a star-shaped domain by the previous proof. See Section 
^, in the proof of Proposition ^.2| for a complementary argument. 
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where the Oi are the holes, that are the bounded connected components of 
lR^\n. A proof can be obtained by analyzing the Ginzburg-Landau equations 
satisfied by a minimizer of G"'^. We finally get : 

MGx,,{u,A)=miG'^-J^{u,A) . (1.12) 



Remark 1.3 . 

If (m, Ae + a) is a solution of the GL-equation then rot a = in the unbounded 



component of J?^ \ri and rot a = const, in each hole (See Lemma 2.1 in | GiPh |). 
It would be interesting to discuss what the possible values of these constants 
are : are they equal to ? 

1.3 Standard results 



The second proposition which is also quite standard (See for example |DGP|) is 



Proposition 1.4 . 

Ifu is a solution of the first GL-equation with the Neumann boundary condition 
then 

\u{x)\ < 1 , \fx e n . (1.13) 

We note for further use that the solutions of the G-L system are in C°°(ri) 
under the assumption that f2 is regular. 



Is the normal state a minimizer ? 



The aim of this section is to give a proof of a result suggested in | BcRu | who 
said " We expect the normal state to be a stable solution for small A...". 
Although, this result is probably known as folk theorem, we think it is useful 
to give a proof (following considerations by M. Dutour in a near context [|Du| ) 
of this property. 



Note th at connected results are obtained in | GiPh | and more recently in [ LuPal |, 
[|LuPa2t . 

Before stating the theorem, let us recall that we have called normal state a pair 
(u. A) of the form : 

iu,A) = {0,A,) , (2.1) 



where A^. is any solution of (|l.l| ). 

We note that this is well defined up to gauge transformation. Moreover, we 
have : 

Lemma 2.1 . 

{0,Ae) is a solution of the GL-system. 

So it is effectively natural to ask if (0, A,,) is a global minimmn. The first result 
in this direction is the following easy proposition about the normal state. But 
let us first introduce : 
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Definition 2.2 . 

We denote by \^^^ the lowest eigenvalue of the Neumann realization in fl of 

-Aa^ := -(V-^A.)' . 
We shall frequently use the assumption 

A'^) > . (2.2) 



Note the following necessary and sufficient condition for this property (cf |He|) 
Proposition 2.3 



The condition (2.£) is satisfied if and only if one of two following conditions is 
satisfied : 

1. is not identically zero in Q,; 

2. is identically zero in 57 but there exists a closed path 7 in 57 such that 

Let us observe that the second case can only occur when J7 is non simply con- 
nected. 

Proposition 2.4 . 

Under condition (2.i.) and if \ g]0,A("'^'[, the pair {0,Ae) is a non-degenerate 



(up to gauge transforms) local minimum of G\^k. 

The Hessian at (0, A^.) of the GL-functional is indeed the map 

5a) ^ ((— A^^ — A) (5m , rot* rot 5a) , 

where we assume that div 5a — Q and 5a ■ v = Q sX the boundary of $7. 
Note that this proof gives also : 

Proposition 2.5 . 

If \ > A*^^^, the pair (0, is not a local minimum of G\,k- 



We refer to [LuPalj for a connected result. Proposition 2.5 does not answer 
completely to the question about global minimizers. The next theorem gives a 
complementary information. 

Theorem 2.6 . 



Under assumption {2_A), then, for any k > 0, there exists Ao(k) > such that, 



for A g]0, Ao(t)], Ga,k has only normal solutions as global minimizers. 



Remark 2.7 



Using a variant of the techniques used in |Du] in a similar context, one can 
actually show that, for any k > 0, there exists Xi{k) > such that, for A G 
]0, Ai(k)], Ga.k has only the normal solutions as solutions of the Ginzburg- 
Landau equations. This will be analyzed in Section |^. 
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Proof of Theorem 2.6: 



Let {u,A) :~ {u,Ae + a) be a minimizer of the (GL) functional. So it is a 
solutionQof (GL) and moreover we have, using ( |1.6[ ), the foUowing property : 



Ga,«(w,A) <0 



(2.3) 



Using the inequaUty — > 



^ and (2.3), we first get, with b — rota : 



— I dx< 



A 



fl2 



(2.4) 



where |n| is the area of fi. 

W e now discuss the hnk between b and a in fl. So we shall only use from 
(0): 



^ J^bUx<^\n\, 

Let us now consider in Q, d the problem of finding a solution of 

rot a ~ b , div d — , 

d ■ I' — , on d^l . 



(2.5) 



(2.6) 



We have the following standard proposition (see Lemma 2.3 in |GiPh|). 
Proposition 2.8 . 

The problem (2.6) admits, for any b G L'^(fl), a unique solution d in H^{fl). 
Moreover, there exists a constant C such that 



,ybeL^ 



(2.7) 



Proof of Proposition |2^ 

Following a suggestion of F. Bcthuel, we look for a solution in the form : 
d — rot*ip. We then solve the Dirichlet problem —Aip = 6 in O. This gives 
a solution with the right regularity. For the uniqueness, we observe that J7 
being connected and simply connected a solution of rot a = is of the form 
d — dd (with 6 G H^{il), and if diva = and d ■ ly on dil, we get the equations 
Ad — and Vd-i^ = on dil, which implies 6 — const, and consequently a = 0.| 



We can now use the Sobolev estimates in order to get 



(2.8) 



^From (2.5), (2.7) and (2.8), we get the existence of a constant C2 such that 

A 



(2.9) 



We actually do not use this property in the proof. 



6 



The second point is to observe, that, for any e g]0, 1[, we have the inequaUty 
|(V - iA)u\^ dx>{l- e)||(V - lAMW^i^n) ~ ^^7^IIHIi^(n) • (2-10) 
Taking ^ — \ and using Holder's inequahty, we get 

|(V - lAM^ dx > ^||(V - - 3||a||i.(,,)||u||i.(,,) . (2.11) 

Using now the eUipticity of — A^^ in the form of the existence of a constant Ci 



L2(n) 



(2.12) 



and again the Sobolev inequahty, we then obtain the existence of a constant C2 
such that 



^ \{V~tA)u\^ dx>(^^~ C2||a||i.(o)) \\iV~iAM\hin)-C'2\\a\\h 



(o)lFllL2(n) ■ 
(2.13) 



We get then from ( ^.Bj ) and (2.9), and for a suitable new constant C (depending 
only on and He), 



A 



(2.14) 



Using the assumption ( |2.2D , this gives w = for A small enough and the proof 
of Theorem 2.6. 

Remark 2.9 . 

Note that with a small improvement of the method, it is possible (taking e — j: 
in (2.10) ) to show that one can choose, in the limit k +00, Ao(k) satisfying : 



Ao(a^) >A«-0(i) 

K 

This will be developped in Section 0. 



(2.15) 



Remark 2.10 . 

Observing that A 1-^ jGx.k.{u, A) is monotonically decreasing, one easily obtains, 
that the set of A's such that (0, A ^) is a global minimum is an interval of the 
form ]0, Aq^'(k)]. Inequality (|2.15| ) implies : 



Ar(«) > A(i) -O(-) 



(2.16) 



Similar arguments are used in [Du| for the Abrikosov's case. We recall that, 
in this case, the domain Q, is replaced by a torus IR^ / C where C is the lattice 
generated over by two independent vectors of IR^ . 

Observing now that k i— > G\^i^{u, A) is monotonically increasing, one easily 
obtains that n 1— > Aq^*(k) is increasing. Using (2.16) and Proposition 3.1, one 
gets that K ^ Aq^*(k) is increasing from to A'-^^for n g]0, +00 [. 
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3 Estimates in the case «; small. 



We have already shown in Proposition 2.5 that, if A > A(^^(k), then the normal 
state is not a minimizer. In other words (see Remark 2.10| ), under condition 
(2.2), we have : 

< A7*(k) < A(i) . (3.1) 

If we come back to the formula ( ]2.14| ) , one immediately obtains the following 
first result : 

Proposition 3.1 . 

There exist constants iiq g]0, A'"'^^] and ao > such that, for A g]0, /ip] satisfying 

A < aoK , (3.2) 
the minimizer is necessarily the normal solution. 

In order to get complementary results, it is also interesting to compute the 
energy of the pair (u, A) = (1,0). This will give, in some asymptotic regime, 
some information about the possibility for the normal solution (or later for a 
bifurcating solution) to correspond to a global minimum of the functional. An 
immediate computation gives : 



dx . (3.3) 

R^ 

We see in particular that when j is small, the normal solution cannot be a 
global minimizer of Ga,k- 

As already observed in Subsection 1.2, what is more relevant is probably the 
integral J^^ dx instead of J ^2 He dx in (3_^). Note also that it would be quite 
interesting to determine the minimizers in the limit k — > 0. We note indeed that 
(1,0) is not a solution of the GL-system, unless Hf, is identically zero in fi. Let 
us show the following proposition. 



-<A-( , (3.4) 



Proposition 3.2 

If 

u- ^ \ . ( 

.^Inll^dx 

and if Q is simply connected then the normal solution is not a global minimum 
Proof. 

Let ij^n a sequence of C°° functions such that 

• < „ < 1; 

• -(/)„ = in a neighborhood of Q; 

• Tpnix) l,Va; ^ 17 ; 
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We observe that 



/ {{1 - ^n)Hef dx ^ [ H^dx. 
We can consequently choose n such that : 



(3.5) 



K<X-( 5 1 , (3.6) 

We now try to find An such that 

• rotA„ = iJjnHe; 

• supp^„ n = 0. 

We have aheady shown how to proceed when is starshaped. In the general 
case, we first cho ose An such that : rot A„ = ipnHe, without the condition of 
support (see ( p^ ) for the argument). 

We now observe that rot A„ = in SI. Using the simple connexity, we can find 
(j)n in C°°{Cl) such that An = V0„. We can now extend outside O as a 
compactly supported C°° function in IB? We then take An = An — V0„. 
It remains to compute the energy of the pair (1, An) (which is strictly negative) 
in order to achieve the proof of the proposition. 

Remark 3.3 . 



In the case when {I is not simply connected. Proposition 3.2 remains true, if we 
replace $7 by f2, where f2 is the smallest simply connected open set containing 

n. 



Remark 3.4 . 

It would be interesting to see how one can use the techniques of |AfDa] for 
analyzing the properties of the zeros of the minimizers, when they are not normal 
solutions. The link between the two papers is given by the relation A = [ndY . 

In conclusion, we have obtained, the following theorem : 
Theorem 3.5 . 



Under condition {2.i), there exists ao > 0, such that : 



1 < < inf Uo, • (3.7) 



.^InH'edx 



4 Localization of pairs with small energy, in the 
case K, large. 

When K is large and A — A^^-* is small enough, we will show as in |Du| that 



all the solutions of non positive energy of the GL-systems are in a suitable 
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neighborhood of {Q,Ae) independent of k > > 0. This suggests that in 
this Umiting regime these solutions of the GL-equations (if there exist and if 
they appear as local minima) wiU furnish global minimizers. Let us show this 
locahzation statement. The proof is quite similar to the proof of Theorem 2.6. 
We recall that we have ( |2.4| )-( 2T0| ). Now we add the condition that, for some 

A<A(i)+r;. (4.1) 

Note that we have already solved the problem when A < A^^-' — ^, so we are 
mainly interested in the A's in an interval of the form [A^^^ — ^, A*^^^ + 77]. 
The second assumption is that we consider only pairs {u,A) e H^{Q)x Hf^^ {IB? ) 
such that 

Ga(w,A)<0. (4.2) 

We improve ( ^.1C| ) into 



2 



(4.3) 

Taking e = we get, using also (|2.9|), the existence of Kq and C such that, for 
A e [0, A'^' + rj] and for n > kq, 

||(V-zA).|P> (^(^(1-^)aW-^)|HP, (4.4) 

for any {u,A) such that G\^k.{u,A) < 0. 



Coming back to (1.1), and, using again the negativity of the energy G\^k.{u, A) 
of the pair (u, A), we get 

\j^\u\Hx<{Ti+^)\\u\\l,^^). (4.5) 

But by Cauchy-Schwarz, we have 

|Mp dx < [ \u\'^ dx)i . (4.6) 

Jn 

So we get 

ML-(n) < C-^)HV + (4.7) 

A hi 

We see that this becomes small with ?/ and It is then also ea sy to co ntro l the 
norm of u in H^{i}). We can indeed use successively ( 2.12| ), ( ^.13| ), ( [4. 2D and 
the trivial inequality: 

||(V-zA)u||i.(f,) <A||u||i.(f,)+GA(u,A). (4.8) 
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The control of {A — A^) in the suitable choice of gauge is also easy through (2.4) 



and (^. 

Note also that if A < A'^^-', we obtain the better 

C 

I|m||l2(o) < — T ■ (4.9) 
kA2 

So we have shown in this section the following theorem: 
Theorem 4.1 . 

There exists ryo > such that, for < rj < rjQ and for A < A^^^ + rj, then there 
exists Kq such that for k > Kq, all the pairs {u,A) with negative energy are in a 
suitable neighborhood 0{ri, i) of the normal solution in H^{Q , C ) x H^{Q, IR^) 
whose size tends to with ri and -. 

Remark 4.2 . 



Using the same techniques as in |Du], one can also show that there are no 



solutions of the Ginzburg-Landau equations outside this neighborhood. This is 
discussed in Section |[ 



5 A priori localization for solutions of Ginzburg- 
Landau equations 



In this section, we give the proof of Remarks 2.7 and 4.S . The proof is adapted 
from Subsection 4.4 in [Du| which analyzes the Abrikosov situation. Similar es- 
timates can also be found in |GiPh] (or in | AfDa]) but in a different asymptotical 
regime. 

We assume that {u, A) is a pair of solutions of the Ginzburg-Landau equa- 
tions (1.3) and rewrite the second Ginzburg-Landau equation, with A = A^ + a 
in the form : 

ia = ^ Im (u • (V - i(Ae + a))u) . (5.1) 



Here L is the operator defined on the space £^^(ri), where, for k G IN* 



E''{n) {a G H'^ifl; iR^) | diva = , a • v/qq = } 



by 



L ~ rot* rot 



-A 



(5.2) 



(5.3) 



One can easily verify that L is an isomorphism from E'^{^l) onto L^(ri). One 
first gets the following 

Lemma 5.1 . 



// {u, Ae + a) is a solution of the GL-system (l.c ) for some A > 0, then we 
have : 

IMI<^. (5.4) 
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Proof of Lemma. 

We start from (5.1) and using Proposition 1.4, we obtain 

iMi'<^ii(v-aHr 

Using the first GL-equation, we obtain : 

A3 



La\\^<\ j \u\\l-\u\')dx 
1^ Jn 



(5.5) 



(5.6) 



Using again Proposition 1.4, we obtain the lemma. 



So Lemma 5.1 shows, together with the properties of L, that there exists a 
constant Cn such that 

\\a\\HHn) < . (5.7) 

This permits to control the size of a when A is small or k is large. In particular, 
using Sobolev's injection Theorem, we get the existence of a constant Cq such 
that : 

< C'q— . (5.8) 

The second step consists in coming back to our solution (u. A) of the Ginzburg- 
Landau equations. Let us rewrite the first one in the form : 



- Aa,u = Am(1 - - 2ia ■ (V - iAe)u - \a\'^u 
Taking the scalar product with u in L'^{fl), we obtain: 



(5.9) 



A II \u\^ ||2 + (-A^^^ , u) < X\\u\\^ + 2\\a\\L^\\u\\^{-AA,u , u) + | |a| li^o 1 P 

<{X+{l + ^)\\a\\l^)\\u\\^+e{-AA^u, u) . 

We have finally obtained, for any e g]0, 1[, and any pair {u, A) solution of the 
GL-equations the foUlowing inequality : 



1 



A / \u{x)\^ dx + {-Aa,u , u) < (A + (1 + 

Jo 1 — e 



1, 



Hp. (5.10) 



Forgetting first the first term of the left hand side in (5.10), we get the following 
alternative : 



• Either u = 0, 

• or 



1-e 



(A+(l + -)||a||i.) 
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If we are in the first case, we obtain immediately (see ( p.l\ ) , the equation La = 
and consequently a = 0. So we have obtained that {u, A) is the normal solution. 

The analysis of the occurence or not of the second case depends on the 
assumptions done in the two remarks, through (5.8) and for a suitable choice of 
e (e = j). So we get immediately the existence of Xi{k) and its estimate when 



K ^ +00. If we now assume (see (|]|)) that A eJA^^) 
back to (5.10) and write : 



V, 



77 [, we come 



A 1^ \u{x)\Ux < • (A + (1 + i)||a||i.) - Ai) \\^ 

Using (O), this leads to 



AIImIM < 



(A + (i + -)IHii.)-AW \n\ 



(5.11) 



This shows, as in (^^), that u is small in with rj and i. 

We can then conclude as in the proof of Theorem 4.1. The control of w in 



is obtained through (5.10) 



Theorem 5.2 . 

There exists 770 > such that, for < 77 < 770 and for A < A*^^^ + 77, then 
there exists kq such that for k > kq, all the pairs (m. A) solutions of the ( G- 
L)-equations are in a suitable neighborhood 0{rj, -i) of the normal solution in 
H\n , (C) X H\n,]R^) whose size tends to with rj and - . 



6 About bifurcations and stability 
6.1 Preliminaries 

Starting from one normal solution, a natural way is to see, if, when increasing A 



from 0, one can bifurcate for a specific value of A. Proposition 2.4 shows that it 
is impossible before A*-^-*. A necessary condition is actually that A becomes an 
eigenvalue of the Neumann realization of — A^_. in il. We shall consider what 
is going on at A^^^ . 

Note here that there is an intrinsic degeneracy to the problem related to the 
existence of an action. We have indeed the trivial lemma 

Lemma 6.1 . 

If (u. A) is a solution, then (exp i6u, A) is a solution. 

This degeneracy is independent of the gauge degeneracy. 
In order to go further, we add the assumption 

A'^' is a simple eigenvalue. (6.1) 

In this case, we denote by ui a corresponding normalized eigenvector. 
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Now, one can try to apply the general bifurcation theory due to Crandall- 
Rabinovitz. Note that, although, the eigenvalue is assumed to be simple, it is 
not exactly a simple eigenvalue in the sense of Crandall-Rabinowitz which are 
working with real spaces. Actually, this is only simple modulo this S'^-action. 
We are not aware of a general theory dealing with this situation in full general- 



ity (see however |GoSc]) but special cases involving Schrodinger operators with 
magnetic field are treated in |0d], | BaPhTa and |Du|. T he a rticle |BaPhTa| 
is devoted to the case of the disk and |Od[ (morc recently [Du|) to the case of 
Abrikosov's states. 

All the considered operators are (relatively to the wave function or order pa- 
rameter) suitable realizations of operators of the type 



u ^ — A^Ti — A/(|up)u , 



with /(O) = 1. 



The main theorem is the following : 

Theorem 6.2 . 

Under the assumptions 



and (6.1), there exist eQ and a bifurcating family of 
solutions {u{-;a),A{-;a),\{a)) in H^{Q.,G)xE^{9)xIR^ , with a G i?(0, eq) C 
(S for the Ginzhurg- Landau equations such that 

u{-; a) = aui + a\a\'^u^^'> a) , with {ui , u^^^) — , 

A{-,a) = Ae + |a|2a2 + |a|4a(4)(.;a) , (6.2) 
A(a) =A(i)+c(K)|a|2 + 0(a4) . 



Here u^-^\-;a) and a^^'^{-;a) are bounded in . 
This solution satisfies, Vs G C , |s| = 1 ; 

u{-;sa) = su{-;a) , A[-;sa) = A{-;ot) 



(6.3) 



Moreover, if c(k) ^ 0, all the solutions (u. A, A) of the Ginzburg-Landau equa- 
tions lying in a sufficiently small neighborhood in x x IR^ of (0, A^, A*-^-*) 
are described by the normal solutions (0,Ae,A) and the bifurcating solutions. 

The constant c(k) will be explicited in the next subsection. 



6.2 About the proof, construction of formal solutions. 

The starting point is the GL-system written in the form 



A(i))ii 



(A - A(i')m - Awjup - 2ia ■ (V - iAe)u - ||a|pw , . 
La = Im (u • (V - iA)u) ^ ' 



We then use the standard method. We look for a solution in the form 

u — aui + ajapua -I- 0(q;^) , 
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a - lap 02 + 0(^4) 

and 

A(a) = A(i'+c(K)|a|2 + 0(a4). 
We can eliminate the iS^-degeneracy by imposing a real (keeping only the par- 



ity). We refer to [Du| for details and just detail the beginning of the formal proof 



which gives the main conditions. We first obtain, using the second equation, 

a2 = -^b2 , (6.5) 

with 

62 := Im (ui • (V - iA^)ui) . (6.6) 
Taking then the scalar product in with ui, in the first equation, we get that 

c(^) = A(i) (^Io - ^Ko^ , (6.7) 

with 



and 



lo := / |ui(x)|^ dx , (6.8) 
Ko = -{ib2-{'V -iAe)ui , ui) . (6.9) 



Remark 6.3 . 

iFrom this expression for c{k), we immediately see that there exists Ki such 



that, for K > Ki, c{k) > 0. Moreover, the uniqueness statement in Theorem 6.1 
is true in a neighborhood which can be chosen independently of k € [ki,-|-cxd[. 

Let us now observe, that, 62 being divergence free, it is immediate by integration 
by part that Kq is real. Computing Re Kq, we immediately obtain : 

Ko^ ReKo^{L-^Ji , Ji) , (6.10) 

where Ji is the current : 

Ji := Im (ui • (V - iAe)ui) . (6.11) 

We observe that Kq > if and only if Ji is not identically 0. In the non 



simply connected case, we shall find a case when Ji = 0. (See Lemma 7.2). 
Following the argument of |Du| (Lemme 3.4.9), let us analyze the consequences 
of Ji = 0. By assumption ui does not vanish identically. If ^1(2:0) ^ 0, 
then we can perform in a sufficiently small ball i?(xo,7'o) centered at xq, the 
following computation in polar coordinates. We write ui — r{x) expi9{x) and 
get Ji = r{x)'^{Ae - V9) = 0. So Ae = V0 in this ball and this implies He = 



in the same ball. Using the properties of the zero set of ui in fl |ElMaQi| and 
the continuity of H^, we then obtain i?e = in il. But we know that, if is 
simply connected, then this implies A^^^ = 0. So we have the following lemma 
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Lemma 6.4 . 

If il is simply connected and X^^^ > 0, then Kq > 0. 

Coining back to the first equation and projecting on uj^, we get : 

U3 = R0V3 , (6.12) 
where V3 is orthogonal to ui and given by : 

V3 := 2a2 ■ {{V - iA,)u,) , (6.13) 
and i?o is the inverse of (—A^^ — A^^^) on the space uj^ and satisfies 

RqUi ~ . 

We emphasize that all this construction is uniform with the parameter P — ^ 
in ]0,/3o]. One can actually extend analytically the equation in order to have a 
well defined problem in [— /3o,/3o]. 

6.3 About the energy along the bifurcating solution. 



The proof is an adaptation of |Du|. Let us just present here the computation of 
the value of the GL-functional along the bifurcating curve. Although it is not 
the proof, this gives the right condition for the stability. For this, we observe 
that if (u, Ag + a) is a solution of the GL-system, then we have: 

GxAu,A) = -^ I + ^ / \mia\^dx. (6.14) 



2 Jn A Jn 

It is then easy to get the main term of the energy of the function for (u, Ag + a) 
with a(-; a) = |Q;pa2(-) + 0{a'^) and u{-; a) — aui{-) + 0{a^). 

GA,«(«(-;a),A(-;a)) - jaj* (^-^ Jju,\^ + ^ JjTota2\^ dx^ + 0(a'^) . 

(6.15) 

Let us first analyze the structure of the term : 

.2 r ^2 



^2 r 

ifi := j |rota2p da; = -^(ia2,a2) . (6.16) 



But we have 



Ki := -^{LbzM) = -^{L^ Ji , Ji) = -^Kq . (6.17) 
With these expressions, we get 
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So we get that the energy becomes negative along the bifurcating solution for 
< I a I < Pq, if the following condition is satisfied : 



(6.19) 



Another way of writing the result is : 
Proposition 6.5 . 



Under conditions {2.i) and Jy , then, if 

JO 

there exists cio > such that, for all a satisfying < |a| < ckq, 
(A(a) - A(i))GA(a),Jw(-; a),A(-- a)) < . 



(6.20) 



(6.21) 



In particular, we have shown, in conjonction with Theorem 4.1, the following 
theorem : 

Theorem 6.6 . 

There exists rj > and kq, such that, for k > kq and A < A^^-* + rj, the global 
minimum of G> is realized by the normal solution for A g]0, A*^^'] and by the 
bifurcating solution for A SjA*^^-*, A*-^-* + rj]. 



In particular, and taking account of Remark 2.10, we have 



Corollary 6.7 . 

There exists Kc such that the map k > Aq^*(k) is an increasing function from 
to A(i) for K e [0, Kc] and is constant and equal to A^^^ for k > Kc- 



Remark 6.8 . 

Note that Theorem 5.2 gives an additional information. For rj small enough and 
K large enough, there are actually no other solutions of the (GL)-equation. 



6.4 Stability 

The last point is to discuss the stability of the bifurcating solution. We expect 
that the bifurcating solution gives a local minimum of the GL-functional for k 
large enough, and more precisely under cond ition ( 6.19). The relevant notion is 
here the notion of strict stability. Following [ BaPhTa , we say that {u, A) (with 
u not identically 0) is strictly stable for Ga.k, if it is a critical point, if its Hessian 
is positive and if its kernel in x is the one dimensional space lR{iu , 0). 
We then have the following theorem : 

Theorem 6.9 . 

Under conditions (2^),(6A), and if ( 6.1!^ is satisfied, then there exists Eq > 0, 
such that, for < |a| < e^, the solution (u(-; a), A(-; a)) is strictly stable. 



We refer to |Du| for the detailed proof. 
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7 Bifurcation from normal solutions: special case 
of non simply connected models. 

7.1 Introduction 

In this section, we revisit the bifurcation problem in the case when is not 
simply connected and when the external field vanishes inside f2. In this very 
particular situation which was considered by J. Berger and J. Rubinstein in 
[ pcRl^ (and later in |HH00l[ , |HH002[ ), it is interesting to make a deeper 
analysis leading for example to the description of the nodal sets of the bifurcat- 
ing solution. The situation is indeed quite different of the results obtained by 
[ ElMaQi[ in a near context (but with a simply connected fl). We mainly follow 
here the presentation in |IIII002| (for which we refer for other results or points 



of view) but emphasize on the link with the previous section. 



7.2 The operator K 



We shall now consider the specific problem introduced by |BeRu] and consider 
the case 

supp i7e n 17 = , (7.1) 
and, in any hole Oi, the flux of satisfies 



(7.2) 



Here we recall that a hole associated to is a bounded connected component 
of the complementary of f2. 



We recall in this c onte xt, what was introduced in [ HHOOl ]. We observe 
that under conditions (7^) and (^^), there exists a multivalued function <j) such 
that expi(j) £ C°°{Ti) and 

d(j) = 2loa , 



where u)a is the 1-form naturally attached to the vector A. 
We also observe that, for the complex conjugation operator F 



Tu = u 



we have the general property 



A_4r 



Combining ( |7.3| ) and (7^), we obtain, for the operator 

K := (exp-i0) r , 

which satisfies 

= Id, 



(7.3) 

(7.4) 

(7.5) 

(7.6) 
(7.7) 
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the following commutation relation 

KAa^AaK. 

Let us also observe that the Neumann condition is respected by K. 
As a corollary, we get 

Lemma 7.1 . 

// V is an eigenvector of , then Kv has the same property. 



(7.; 



This shows that one can always choose an orthonormal basis of eigenvectors Uj 
such that Kuj — Uj. 

7.3 Bifurcation inside special classes. 



Following |BcRu| (but inside our point of view), we look for solution of the GL 

(7.9) 



equation in the form (u, A^) with Ku = u. Let us observe that 
Lj^ifl; €) ■.= {ue L^{n; (C ) | ifu = u} 



is a real Hilbert subspace of L'^{^1; (U ). 

We denote by H'^ the corresponding Sobolev spaces 



(7.10) 



We now observe the 



Lemma 7.2 . 

If u ^ Hj^, then Im (fi • (V — i^e)w) — almost everywhere. 



Proof of Lemma 7.2 



Let us consider a point where u 7^ 0. Then we have u — pexpiO with 29 = cj) 
modulo 2-kZ. Remembering that Ae = ^V^, it is easy to get the property. | 

Once this lemma is proved, one immediately sees that {u, A,,) (with Ku = u) is 
a solution of the GL system if and only if u e H]^ and 



-Aa^u — Au(l 



, 



[V -iA^)u-v , on 00 



(7.11) 



We shall call this new system the red uced G L-equation. But now we can apply 
the theorem by Crandall-Rabinowitz |CrRa | . By assumption (O), the kernel of 
(— Aa„ — A^^^) is now a one-dimensional real subspace in L"]^. Let us denote by 
ui a normalized "real" eigenvector. Note that ui is unique up to ±1. Therefore, 
we have the 

Theorem 7.3 . 



Under assumptions (6.1), (7.1) and (l.i.), there exists a bifurcating family of so- 
lutions 
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a), A(a)) in H]^ x IR^ with a £] — eo,+eo[, for the reduced GL-equation 
such that 

u{a) = aui + v{a) , 
(mi , v{a))L2 ^ , 
\\v{a)\\H2^n)=Oil) , 



with 



Moreover 



A(a) = A(i) + ca^ + ©(a^) 
c = A« 



|ui|^ dx 



u(— a) = —u{a) , A(— a) = A(a) 



(7.12) 
(7.13) 
(7.14) 

(7.15) 



Remark 7.4 . 

Note that the property (7.15) is what remains of the S'^-invariance when one 
considers only "real" solutions. 

Let us give here the formal computations of the main terms. If we denote by 
Lq the operator La := — A^^ — A^^^. Writing v{a) = 1*3 + 0{a), we get : 



Projecting on m, we get (7.14). Projecting on uj^ and denoting by i?o the 
operator equal to the inverse of Lq on this subspace and to on Ker Lq, we get 



U3 = -A«i?o(ui|"in = -A(i)i?o(wi|wip-cui) 



(7.16) 



Remark 7.5 . 

By the uniqueness part in Theorem 6^, we see that the solution {u{-; a), Ae) is 
actually the solution given in this theorem. 



Another remark is that 



aw 

2 



a\ \ui{x)\Ux) + Oi\af) , (7.17) 
Jn 



so that when a 7^ the energy is decreasing. This is of course to compare with 
(3.16) (note that we have = 0). Once we have observed this last property, 
the local stability of the bifurcated solution near the bifurcation is clear. 

The second result we would like to mention concerns the nodal sets. In the 
case when fl is simply connected, the analysis of the nodal set of u when (w. A) 
is a minimizer of the GL-functional is done in [ElMaQij, using the analyticity 
of the solutions of the GL-equation and techniques of Courant. 

In the non si mply c o nnected c ase, very few results are known. The following 
theorem is true 1BcRu[ , lHH002i : 
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Theorem 7.6 . 

Under assumptions (2.i ), (7.1) and (7.i), there exists ei > such that, for any 
a €]0,ei], the nodal set of u{a) in Hj^ slits Cl in the sense of [HHOOlj. In 
particular, if there is only one hole, then the nodal set of u{a) consists exactly 
in one line joining the interior boundary and the exterior boundary. 



An elegant way to recover these results (see [HH001|, |HH002|) is to lift 
the situation to a suitable two- fold covering fi^. 
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A Analysis of the various scalings. 



When considering asymptotical regimes, it is perhaps useful to have an inter- 
pretation in terms of the initial variables. According to the statistical inter- 



pretation of the Ginzburg-Landau functional (See for example [ BeRuSc |), the 
starting point is the functional {v,A) J-{v,A) with : 

^ Iroti - i/eP dx 
+ L£\i^-^'-FA)u\'dx 
+ J^^{a\u\^ + I\u\^) di. 

Here a is a parameter which is proportional to (T — Tc) (we are only interested in 
the case a < 0) and b is essentially independent of the temperature. The other 
parameters are standard : h — h is the Planck constant, e is the charge of 
the electron and m is the mass of the electron. With u = tKu and A — —A, 

\a\ c ' 

we obtain : 

with He — j^He, ^ — and k = ^(^)^- Here we emphasize that between 

the two functionals, no change of space variables is involved. 

Let now compare with another standard representation of the Ginzburg- 
Landau functional. We make this time the change of variables x — and 
if we change u and the 1— form corresponding to A accordingly, we obtain the 
standard functional : 

£{u,A) = GxA'^^A) , 

with 

E{u,A) = j^^{-\u\^ + \\u\'^)dx 
+ /o l(V - iA)u\'^dx 
+ |roti - Hel'^dx , 



with 



He — ^He 



Here we observe that the open set il is not conserved in the transformation. We 
have to keep this in mind when comparing in the limit k — > -t-oo the contribu- 



tions of Sandier and Serfaty |SaSe| or |LuPal with the results presented in this 
paper. 
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